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MODELING  MAINTENANCE  DATA  BY  EXPONENTIAL,  WEIBULL,  GAMMA 
NORMAL  AND  LOGNORMAL  DISTRIBUTIONS 

1 . INTRODUCTION 

1.1  Objective 

Messrs.  Blanchard  and  Lowery's  1969  study,  listed  as  reference 
number  two.  Indicated  that  maintenance  data  model  to  the  normal  and 
lognormal  distributions  as  opposed  to  some  of  the  other  frequently  used 
distributions.  This  report  evaluates  a sample  of  Army  vehicle 
maintenance  data  to  determine  how  well  It  satisfies  their 
conclusion. 

1.2  Background 

The  time-to-repair  data  used  In  this  analysis  were  derived  from 
maintenance  data  which  were  taken  on  an  Army  M151A1  1/4  ton  utility 
truck.  These  data  were  collected  through  the  Army  Integrated  Equipment 
Record  Maintenance  Management  System  (TAERS).  Several  hundred  time-to- 
repair  unscheduled  corrective  actions  were  recorded  from  an  Army 
vehicle.  Of  these,  a random  sample  of  366  samples  was  utilized  for 
this  evaluation.  Isolation,  removal,  replacement,  alignment  and 
verification  times  are  not  considered  In  this  data  sample  analysis. 

The  personnel  responsible  for  restoring  a piece  of  equipment  to  Its 
original  operable  state  were  Individuals  with  some  formal  training  to 
be  mechanics.  History  of  these  mechanics  who  record  time-to-repair 
data  Indicated  that  they  tended  to  round  off  their  restoring  times  to 
the  nearest  one  tenth  hour.  This  will  be  discussed  In  more  detail  In 
a later  section. 

Maintainability  time  can  be  defined  as  "a  characteristic  of  design 
and  Installation  which  Is  expressed  as  the  ability  to  retain  or  restore 
an  item  to  a specified  condition  within  a given  period  of  time." 
Ordinarily,  these  times  are  considered  to  be  continuous  random  variables; 
however,  these  data  appeared  to  be  recorded  by  rounding  off  to  the  near- 
est one-tenth  hour  and  possibly  rounded  to  the  nearest  one-half  hour. 

As  a result,  raw  data  In  Table  I (Page  28)  based  on  random  sample  times 
will  be  analyzed  using  discrete  random  variables  (x)  since  for  each 
outcome  (x^),  there  Is  a P(x^)  ^0,  1 = 1,  2 ...  n. 

Figure  I Is  a histogram  of  corrective  repair  times  in  hours  based 
on  the  non-grouped  data  In  Table  II.  Figure  II  Is  an  Identical  histo- 
gram to  Figure  I except  Its  Interval  widths  were  Increased  to  adjust 
for  round  off  In  recording  repair  times.  These  round-off  data  correspond 
to  adjusted  Intervals  In  Table  III. 

The  fitting  of  these  data  to  the  widely  used  exponential,  Welbull, 
ganma,  normal,  and  lognormal  distributions  are  analyzed  In  respective 
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sections  of  this  report.  Although,  this  report  will  show  computations 
for  both  non-grouped  data  (Table  II)  and  grouped  data  (Table  III),  the 
conclusion  will  be  determined  from  adjusted  interval  data  given  in  Table 
III.  Each  section  utilizes  the  same  basic  interval  widths,  even  though 
the  forms  of  the  exponential,  Weibull,  gamma,  normal,  and  lognormal 
distributions  are  different. 

A statistical  comparison  test  using  observed  data  with  theoretical 
results  of  the  assumed  distributions  had  to  be  selected.  Such  test 
is  called  the  "Goodness-of-fit  test;"  two  are  used  in  this  analysis. 

2 

These  tests  are  the  Chi-square  (x  ) and  Kolmogorov-Smirnov  (K-S)  tests. 

2 

The  X goodness-of-fit  test  can  be  used  with  large  samples  to  test 
the  validity  of  any  assumed  discrete  or  continuous  distribution.  The 
Kolmogorov-Smirnov  is  useful  in  large  or  small  samples,  however,  it  is 
restricted  to  testing  the  continuous  distributions.  Both  goodness-of-fit 
tests  were  applied  in  this  analysis,  although  the  Kolmogorov-Smirnov  (K-S) 
test  analysis  is  not  shown  in  this  report.  Figure  III  gives  a pictorial 
representation  of  additional  reasoning  for  not  choosing  the  K-S  Test. 

It  is  easily  seen  in  this  figure  that  if  a curve  were  fitted  through 
the  center  of  maintenance  points  in  the  0.5  interval,  the  standard 
deviation  would  be  at  least  0.2  as  compared  to  a smaller  number  at  a 

2 

significant  level  of  0.10,  0.05,  or  0.01.  Consequently,  the  x 
goodness-of-fit  test  was  decided  to  be  the  better  of  the  two  tests. 

II.  EXPONENTIAL  DISTRIBUTION 

Table  I presents  366  individual  corrective-maintenance 
repair  times  for  all  kinds  of  truck  failures.  Such  data  will  be  used 
in  the  assumed  underlying  distribution  to  compute  an  estimate  of  its 
parameter.  The  maximum  likelihood  method  is  used  in  the  following 
equation  to  estimate  the  e parameter  of  an  exponential  distribution 
(l-EXP(-X/0)) 


e 


e 

where 

6=  Mean  time  to  repair 
N = Sample  size 

X = Maintenance  repair  time  (Hours) 
^ = Failure  rate 


= I X.  = - = Sample  Mean 


= 242.8/366  = .6634 
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With  the  estimated  parameter,  a statistical  goodness-of-fit  test  is 
used  to  evaluate  the  assumed  distribution.  Such  a goodness-of-fit  test 

2 

is  the  Chi-square  (x  ) test.  Three  important  observations  about  the  use 
of  this  test  are  as  follows: 


1.  It  may  be  used  with  a discrete  or  continuous  distribution. 

2.  It  allows  the  use  of  estimation  of  parameters  and  underlying 
distribution. 


3.  It  requires  "large"  sample  size  (i.e.,  20  or  more). 


The  steps  necessary  to  apply  the  Chi-square  goodness-of-fit  test  are 
given  below: 

1.  Partition  the  range  of  the  variable  into  intervals  with  each 
interval  closed  from  the  left  side. 


2.  Determine  the  number  of  sample  observations  falling  within  each 
interval . 

3.  Determine  the  level  of  significance,  which  is  defined  as  the 
risk  of  rejecting  the  underlying  distribution  if  it  is,  in  fact,  the 
real  distribution. 


Compute 


i=l 


'i 


where 

0^.  = Number  of  sample  observations  in  the  i^*^  interval 
= Expected  number  of  observations  in  the  i^*’  interval 
N = Number  of  intervals 


2 2 

5.  Determine  if  x is  greater  than  N P 1 
where 

P is  number  of  parameters  estimated  from  the  data. 
N is  number  of  interval  lengths 
o is  risk  level 
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2 2 2 ? 

^o,  N-P-1  found  in  a x table.  If  x is  larger  than  x^  N-p-i 

a decision  is  made  to  reject  the  distribution  under  test,  otherwise,  we 
do  not  have  sufficient  evidence  to  reject  the  assumed  distribution. 

The  use  of  original  repair  times  data  without  adjusting  interval 
length  is  unrealistic  since  the  recorders  apparently  rounded  off  these 
repair  times.  An  analysis  of  original  data  with  one-tenth  hour  intervals 
(Table  II)  was  used  to  verify  this  conjecture.  Subsequently,  some  of 
these  interval  lengths  (Table  II)  were  extended  to  include  more  sample 
observations  in  each  of  these  intervals  (Table  III).  The  data  show 
that  the  recorders  tended  to  round  the  actual  time  it  takes  to  repair 
an  item  to  the  nearest  12  minutes  (.2  hour),  30  minutes  (.50  hour),  and 
1 hour.  Table  II  shows  that  out  of  366  corrective  repair  maintenance 
times,  there  are  54,  148  and  67  of  such  occurring  at  .2,  .5  and  1.0  hour, 
respectively.  The  selection  of  the  interval  widths  for  the  Chi-square 
test  was  based  on  the  rounding  phenomena  since  the  expected  number  of 
observations  in  any  interval  must  be  at  least  five  for  successful  use 
of  the  Chi-square  test,  width  2.25  to  4.05  was  selected.  The  risk  level 
(a)  of  .05  using  Chi-square  test  is  generally  used  in  rejecting  the 

assumed  distribution  if  the  assumption  is  true.  The  x value  for  this 
distribution  is 

2 _ 2 

N-P-1  ■ ^.05,  18 

= 29 

where 

a = .05 
N = 20 


P = 1 

2 

based  on  unadjusted  data  interval  used  in  Table  II.  Also,  the  x value 
of 

^a,  N-P-1  = x^05,  5 
= 11.07 

where 

a = .05 
N = 7 
P = 1 


I 

i 
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is  based  on  the  grouped  data  intervals  in  Table  III.  The  goodness-of-fit 
test  uses  these  tables  Chi-square  values  of  29  and  11.07  to  compare  with 

N (o.-ej2  N (O.-E.)^ 

the  calculated  t \ ^ of  977  (Table  IV)  and  F p - of  73  Table 
i=l  “"i  2 ^ 

V,  respectively.  As  a result,  both  values  are  less  than  both  table 
values  which  indicates  that  the  assumed  distribution  is  not  exponentially 
distributed. 

III.  WEIBULL  DISTRIBUTION 


The  Weibull  distribution,  considered  in  this  section,  is  a 
possible  model  for  repair  times.  Such  a distribution  is  given  by  the 
following  equations: 


F(X)  = 


1-EXP  {-[h 


X>0 

n>0 


0 otherwise 


where 

X = Maintenance  repair  time  widths 
n = The  scale  parameter 
6 = The  shape  parameter 

Some  of  the  same  procedures  are  applied  to  this  distribution  as  the 
exponential  distribution  as  noted  above.  Applying  the  method  of  match- 
ing moments  to  estimate  the  parameters,  n and  b are  calculated  from  the 
following  respective  equations  using  previously  calculated  estimations  for 
mean  e and  standard  deviation  o. 

A 

^ 0 

n = — T is  the  scale  parameter 

Hfnr 

where 

6 = 1/b  is  the  shape  parameter 
or  b = 1/6 

The  b is  calculated  from  equation 
2 

aO  a 

1 + £ = liLJjJ 

i 2 


n 


1 + £ = 2r(2b) 

b (r(b))^ 

1 + (.6021)^  = 2r(2b)) 
(.6634)^  b(r(b))^ 


1.82  = 

b(r(b))2 


which  must  be  solved  first  by  iteration  below  and  using  a gamma  function 
table. 


b 


2r(2b)) 


b(r(b))^ 


2r(2b) 

b(r(b))2 


1 


.9 


2r(2)  = 2 

2r(1.8)  = 
1.8628 


l(r(l))^=l 

.9(r(.9))2= 
1 .0278 


2 

1.82 


Therefore  b = .9  satisfies  the  above  relationship,  now  an  estimate  for 
6 is  calculated  below. 


i 

.9 


6 = 1.11 


and  the  estimate  for 

- _ .6634 
^ ■ r{1.9) 

_ .6634 
' .96l>7 


= .6898 

The  other  change  is  the  computation  of  the  expected  number  of  frequencies 

e 

E.  = N (e  e ( ; ) ) 

where 

6 and  n defined  above 
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X^.  = Intervals  widths 

N = Number  of  observations  = 366 
2 

The  X values  for  this  distribution  differ  from  the  exponential 
distribution.  The  two  parameter  Weibull  distribution  is  used  to 

calculate  ^-P-l  ^05. 17  = 

and  X Q5  ^ = 9.5  given  in  respective  data  tables  II  and  III.  Note  that 
these  values  for 


N (0.-E.)‘ 
ZQ  < I ’ ^ 


i=l  ®i 


876 


and 


N (0.-E.)‘ 
9.5  < I — ^ 
i=l  ''i 


54 


imply  that  the  assumed  distribution  is  not  a Weibull  distribution. 


IV.  GAMMA  DISTRIBUTION 


The  gamma  distribution  is  applied  to  these  data  to  test  how  well  such 
data  fit  this  assumption.  This  analysis  is  similar  to  those  presented 
earlier.  A method  for  estimating  the  parameters  6 (scale  parameter)  and 
a (shape  parameter)  is  the  maximum  likelihood  procedure  based  on  the 
gamma  density  function. 

(1)  f(X;3.  a)  = (3“  X“"''  e‘^^)/r(a) 

the  natural  log  of  t is  L 

N 

where  t = ir  f(X*,3,a)  a is  the  shape  parameter 

i=l  3 is  the  scale  parameter 

N 

(2)  L = a N In  3 - N In  (r(a))  + (a-1)  1 In  (X.) 

i=l  ^ 

N 

- e I X^ 
i=l  ^ 

(3)  N In  3 - N (a)  + In  (X^)  = 0 
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(4) 


ik  - oN 

36  " B 


M 

I X. 

i=1  ^ 


Setting  the  partial  derivatives  in  equations,  2 and  4 equal  to  0. 
where  (|  (a)  is  In  (r(ci)) 

equations  (3)  and  (4)  can  only  be  solved  by  trial  and  error,  therefore 
equation  (4)  becomes 

. fi  N 

(5)  ct  = f I X.  = 6 X 
^ i=l  1 

where 

N = Total  Number  of  Samples 
X^.  = Individual  Sample 
X = Mean  of  Sample  = .6634 
6 = Scale  Parameter  = 4.32186 
a = Shape  Parameter 
Therefore  a is  estimated  to  be 
a = 4.32186  (.6634)  = 2.86712 
and 

1 N 

(5)  6 = exp  {^{a)  - ^ ^ In  (x.)) 

^ i=l  ^ 

which  is  simplified  to 
j . k e 

where  N 

. i I In  (x.) 
k = e N i=l  ^ 

The  approximation  for  this  is 


;<^(ex)  , (g  X - i)  (1  + — J ) 


j(^(6x)  ^ (3X.1/2)  + 


24(eX-l/2)' 
1 


24(bX-1/2) 
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Therefore,  the  equation  for  e in  (5)  above  resolves  to  a quadratic  equation 

ae^  + be  + c = 0 

Hence,  the  positive  root  of  the  above  equation  is 


2a 


Therefore 

B = 4.32186 


These  parameters  may  also  be  obtained  by  the  method  of  matching  moments, 
however  it  does  not  yield  the  minimum  variance  estimate  for  large  samples, 
therefore,  the  maximum  likelihood  method  was  utilized.  As  a result, 
these  estimated  parameters  and  repair  times  are  used  in  the  underlying 
distribution  analysis.  Although  the  area  under  this  distribution  curve 
is  not  computable  due  to  its  non-closed  form,  it  can  be  approximated  using 
Simpson's  Rule  applied  to  the  gamma  density  function  (g.d.f.)  that  is 


f(x) 


r(ct) 


X>0 


Tables  VIII  value  of  605  and  Table  IX  value  of  80  give  results  of 
computerized  g.d.f.  calculations  using  original  data  intervals  and 


adjusted  data  intervals  respectively.  The  computer  values  of 


I (Oi-^i)' 


^i 


i=l 

2 2 

in  both  tables  above  are  compared  To  x and  x gg  3 for  values  of 

605  and  80,  respectively  It  is  concluded  that  the  Chi-square  test  result 
shows  that  this  is  not  a gamma  distribution. 


V.  TRUNCATED  NORMAL  DISTRIBUTION 


The  normal  distribution  is  the  best  known  distribution  of  statistics. 
It  would  certainly  be  appropriate  to  investigate  maintenance  times  with 
the  aid  of  the  normal  distribution.  Unfortunately,  this  distribution 
has  a characteristic  which  at  times  precludes  its  use  with  certain  kinds 
of  data.  Essentially  the  normal  distribution  is  defined  for  a random 
variable  spanning  the  whole  real  line.  The  maintenance  times  now  under 
study  can  never  be  negative  and  are  best  described  by  a random  variable 
whose  range  is  the  positive  half  line.  Still  it  is  possible  to  use  the 
normal  distribution  in  the  investigation  of  maintenance  times  if  the 
forbidden  region  of  the  maintenance  time  random  variable  has  a low 
probability  of  occurrence  when  described  as  the  appropriate  normal 
distribution.  A rough  rule  of  thumb  to  judge  the  applicability  of  the 
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normal  distribution  to  positive  data  is  to  use  this  distribution  only 
if  the  ratio  of  the  sample  mean  to  the  sample  standard  deviation  exceeds 
three.  In  our  case  this  criterion  is  not  satisfied.  The  mean  of  the 
sample  data  is  .66  and  the  standard  deviation  estimate  is  .60  based  on 
the  following  equation: 

n X. 

0 = y -—sanole  mean 
i=l  " 

242.8 
" 366 

= .6634 


where  x^  = each  corrective  repair  time 
n = Total  number  of  observations 


n = 366 


where  x^.  and  n = same  as  above, 


therefore  o = 


8432.56 

133590 


o = .6021 


A normal  variate  with  these  values  for  mean  and  standard  deviation  would 
have  a probability  .14  of  being  negative.  This  is  much  too  high  a prob- 
ability to  be  associated  with  data  values  which  cannot  possibly  occur. 

In  order  to  circumvent  this  difficulty  the  normal  distribution  may  be 
truncated.  In  effect  the  area  under  the  normal  density  curve  to  the  left 
of  some  point,  the  truncation  point,  is  ignored.  The  remaining  part  of 
the  density  function  is  multiplied  by  an  appropriate  constant  to  insure 
the  total  area  under  the  curve  is  one.  The  truncated  normal  distribution 
is  specified  by  either  two  or  three  parameters,  depending  on  whether  the 
truncation  point  is  known.  In  the  present  instance  the  variable  of  interest 
is  bounded  by  zero  and  thus  there  are  only  two  parameters.  If 

w = Mean  Estimate  and 

a = Standard  Deviation  Estimate  of  the  untruncated  normal  distribution, 
then  the  density  function  f(x),  of  the  truncated  distribution  is 
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f(x)= 


where 


and 


— - „-l/2 


G is  defined  by 
_L  Z _ 


1 


1 


' - '^(Z)  = '■  I- 


ds 

2 


$2  = 


x>0 

x^O 


G(Z)  is  the  cumulative  probability  function  for  the  standard  normal 
variate.  It  is  convenient  to  define  as  an  auxiliary  parameter  the 
standardized  point  of  truncation  TS.  That  is 


TS  = (^) 

x = 0 for  this  case 
therefore  T = - 

o 

It  should  be  clearly  understood  that  even  though  three  parameters 
will  be  estimated  only  two  are  needed  to  specify  the  distribution. 

The  method  of  matching  moment  will  be  used  to  estimate  the  parameters. 
According  to  Cohen,  this  procedure  is  equivalent  to  the  maximum  likelihood 
method,  and  in  addition  is  easily  described.  Furthermore,  along  the  way, 
expressions  for  the  population  mean  and  population  standard  deviation  will 
be  derived.  The  derivation  begins  with  the  computation  of  the  first 
moment  (M-j)  and  second  moment  (M2)  of  the  distribution. 

/ se-^/2  (— ) ^,5 

TS 


TS 


«,  . E(x)  = J 


1 


^ a (l-G(TS)) 


M2  = E(x^)  * ^ ^ 

^ ^ a (l-G(TS)) 
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To  simplify  : 


Let  V = ^ 

0 

d,  = * 


- * 

j (aV+u)  ® ^ dv 
TS 


By  integration 

M-i  = — ^ T-^ 

As  above  let 


a e-  nsiS 


1 Z .2 

S(Z)«^  /_e-|  ds 

and  define  L(x)  by 

x2 

' L1-G(x)J  ‘ ® ^ 

By  substitution  of  L(x) 

= a L(TS)  + u = Population  Mean 


J ^ 

[l-G(TS):  ^ 


Now  simplify  M2: 


M = j_  1 r o . v2 

2 D-G(TS)  3 (y  + aV)  e 7 


( V 

+ ZCM  v+o^v^)  e"  2 


By  substitution  of  G(TS) 


l-G(TSn  . 2(au)^S  ''®  ^ dv  . 4s  e"  7 dv 


^[1-G(TS)1  /27:1-G(TS)] 


By  integration 


2 . iiy.  2 2 / -t 

M2  = ji  + 2 e 2 ^ g TS  e 2 ^ a TS  e 7 

^ :i-G(TS)]  Cl-G(TS)]  [l-G(TS)] 

Again  by  substitution  of  L(TS)  and  G(TS) 

+ 2 ayL(TS)  + / (TS)  L(TS)  + 


To  obtain  the  estimate  for  y,  a and  TS  the  population  first  and  second 
moments  are  replaced  by  the  corresponding  sample  values.  For  this  purpose 
the  equations  for  y,  a,  and  TS  may  be  written  as 

1^2  - = y^+2ay  L(TS)+o^(TS)L(TS)+a^  - (oL(TS)  + y)^ 

M2  - M^^  = y^+2ay  L(TS)L(TS)+a^  - o^L^(TS)  - 2ay  L(TS)  - y^ 


= [g  L(TS)  + yf 
M^^  = g^  L^(TS)  + 2gyL(TS)  + v 


Note  that:  y = - oTS  = -1.7867  (2.0173)  = -3.6043 

Where  y = associated  normal  mean  and  a associated  Normal  Standard 

Deviation 


M^^  = g^  D.^(TS)  - 2 TS  L(TS)  + (TS)^] 


* ii2 

”1 

(TS)  - 2 TS  L (TS)  + (TS)^ 
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1. 

Since  L(x)  is  not  in  closed  form  and  an  intergral  part  of  equations 

A A 

and  Mg,  a computer  program  is  used  to  calculate  estimates 

A A A A A A 

for  Mp  Mg,  and  for  y and  a,although  u and  a are  completed  later.  Now 

the  ratio  estimate  for 


.3615  _ 
.4400996  " 


.8214 


Where  M^  = .6634 


Mg  = .80164 


is  computed. 


These  two  parameters  are  now  used  in  the  truncated  normal  probability 
density  function  (P.d.f.)  below  to  estimate  the  area  under  this  assumed 
distribution  curve. 


P.d.f.  f(x) 


1 1 


n-G(x)D 


-1/2 


2^ 


Since  this  distribution  is  not  in  closed  form,  it  requires  the  use  of  a 
standard  normal  table  to  compute  the  area  under  the  truncated  normal 
distribution  curve.  Calculated  results  of  such  areas  are  generated  in 

N (0.  - E.)2 

the  computation  of  ) f — ■ — in  Tables  X and  XI  for  unadjusted 

i=l  H 

intervals  and  adjusted  intervals,  respectively.  From  Table  X we  obtain 
N (0.  - E.)2 

I E ® which  is  greater  than  x Q5  ig  = 26  and  from 

N (0.  - E.)^  2 

Table  XI  I ? = 68  which  is  greater  than  x nt;  a ~ 9.5. 


The  goodness  of  fit  test  indicates  that  the  observed  maintenance 
times  do  not  come  from  a normal  distribution. 


VI.  LOGNORMAL 

The  lognormal  distribution  requires  that  data  be  analyzed  in  similar 
manner  as  in  previous  sections.  However,  the  natural  log  of  individual 
maintenance  repair  times  data  from  Table  I is  used  to  calculate  the  two 

parameters  which  are  the  sample  mean  (e-ip)  and  standard  deviation  (a-i^)* 

1.  ©Ip  * ^ ^~N^^  ” Sample  mean 
where 

N = Total  number  of  Maintenance  Repair  Times 
N = 366 

X^  = Individual  Maintenance  Repair  Time 

then 
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" _ -257.55768 

®ln  “ 366 


0^^  = -.7037 


and  1 5 

« ^/Nldn  - (^In  X.) 

°ln  1 NTrTTj 


= Sample  Standard  Deviation 


where 

N(N-l)  = 366  (365)  = 133590 
Nl  (In  X^)^  = 145614.07 

(Iln  X^)^  = 66335.958 
yields 

-In  - .7703 

The  lognormal  distribution  function  is  also  not  in  closed  form,  however, 
its  area  under  this  curve  can  be  approximated  by  similar  methods  as  were 
used  for  the  normal  distribution  approximatio.i.  First,  a change  of 
variable  is  made,  that  is  from  lognormal  density  function  (l.d.f.). 


f(x) 


1 


-1/2  (- 


In  X- 


-In 


^In 


XaV^rT 
let  y = In  X 

- i 

clx  ~ X 

Then 

dx  _ 

37-  " 

f(x)  transforms  to  g(y)  = f(x)l^l 

' 2 

1 -1/2 

g(y)  “ e ' ' " Ms  standard  normal  density  function 

which  is  identical  to  normal  density  form. 
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• . * .. 


Now,  the  area  under  the  lognormal  distribution  can  be  found  using  the 
normal  distribution.  Finally,  the  Chi-square  test  at  the  same  .05 
risk  level  is  used  to  determine  if  the  assumed  distribution  is  lognormal. 

N (0.  - 

Tables  XII  and  XIII  show  computerized  results  of  > f — ■ — 

i=l  ''i 

based  on  the  same  data  intervals  used  in  previous  section.  Table  XII 
18  (0.  - EJ2 

value  of  706  for  > — —f — = — and  Table  XIII  value  of  12.79  are  compared 
i=l  “"i 
2 2 

’^.OS,  15  ^.05,  4*  '"espectively. 

These  comparisons  contributed  the  results  of 
2 

X =25  which  is  less  than  706 

.05,  15 

2 

X 05  4 = 9.49  which  is  less  than  12.79  which  infer  that  such 

lognormal  distribution  assumption  is  invalid  at  the  .05  risk  level. 
However,  for  this  particular  distribution  another  risk  level  (.01)  can 
be  used  to  satisfy  the  Chi-square  test  based  on  the  adjusted  intervals 
due  to  rationale  stated  in  the  introduction.  Such  information  is  shown 
2 

by  X 01  4 = 13.28  which  is  greater  than  12.79.  Therefore,  this  means 

that  at  the  .01  risk  level  there  is  not  sufficient  evidence  to  say  that 
this  distribution  is  not  lognormally  distributed. 


VII.  CONCLUSIONS 

The  analytical  results  of  each  section  are  compared  to  determine 
if  the  discussed  distribution  can  be  used  to  model  these  data.  Table 
XIV  gives  a summary  of  calculated  data  ratios  which  are  compared  to  the 
2 

X goodness-of-fit  test  at  the  5 percent  risk  level.  It  is  seen  that 
each  distribution's  computed  ratio  yielded  a larger  value  than  its 

O 

corresponding  x value  for  unadjusted  as  well  as  adjusted  intervals. 
However,  the  lo-^normal  distribution  does  come  closest  to  satisfying 
such  tests.  T'»refore,  based  on  this  analysis  the  lognormal  distribution 
is  the  best  of  these  five  models  for  the  maintenance  data  analyzed. 
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TABLE  I 


1 

Maintenance  Repair 

Times  (Hours) 

.1 

.5 

.2 

1.0 

.5 

.5 

.5 

3.0 

.5 

.2 

4.0 

.1 

.1 

.2 

.5 

.2 

.5 

.5 

.5 

1.0 

.5 

.5 

.3 

1.0 

1.0 

.2 

.2 

.5 

.2 

1.0 

.5 

.5 

.5 

.5 

.2 

.5 

.5 

1.0 

.2 

.3 

.5 

.3 

1.0 

1.0 

.5 

.5 

1.0 

.2 

.5 

.1 

l.o 

.5 

.3 

.3 

.2 

.2 

.4 

.5 

.5 

.2 

1.0 

.2 

.5 

.7 

.1 

.5 

.2 

.1 

.5 

.5 

.5 

.5 

.5 

.2 

.5 

.5 

.6 

.1 

.5 

.5 

.4 

.5 

1.0 

.5 

.5 

.5 

.2 

.6 

1.0 

.5 

.2 

.2 

4.0 

.6 

.5 

.5 

.5 

.5 

2.0 

1.5 

.6 

.5 

1.0 

.2 

.1 

.1 

.1 

1.0 

.5 

.5 

.5 

1.0 

.5 

.2 

.5 

.6 

.3 

3.5 

.1 

.2 

1.5 

.5 

.5 

.5 

1.0 

.5 

.3 

.3 

2.0 

.4 

2.0 

.2 

.3 

.5 

.5 

1.0 

.5 

4.0 

.5 

.2 

.4 

.5 

1.3 

.2 

.5 

.5 

1.0 

.5 

.5 

.5 

1.0 

.5 

.5 

.6 

.4 

.5 

.2 

.5 

1.0 

.5 

.5 

1.0 

.5 

2.5 

1.0 

.2 

.5 

1.5 

.3 

.5 

.5 

.5 

.5 

.5 

.5 

.5 

.5 

.5 

.3 

.5 

1.0 

.3 

.5 

.1 

1.0 

.5 

1.0 

.5 

.5 

.5 

.2 

.2 

2.5 

2.0 

.7 

.5 

2.0 

1.0 

.5 

1.0 

.5 

1.0 

.2 

.5 

.5 

1.0 

.5 

.1 

1.0 

.5 

.1 

.5 

.5 

.5 

1.0 

.5 

.5 

.2 

.7 

1.0 
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TABLE  I (Continued) 


Maintenance  Repair  Times  (Hours) 


.1 

.5 

.5 

.5 

.5 

4.0 

.5 

,2 

1.5 

.5 

.2 

1.0 
1.0 

.2 

.5 


1.0 

2.0 

1.0 

.5 

.5 

.5 

.8 

.9 

.2 

1.0 

1.0 

.1 

.6 

1.0 

.2 


.2 

1.0 

.2 

2.0 

1.5 

.5 

.5 

.5 

.5 

.5 

.5 

.6 

.5 

1.0 

1.0 


1.5 

.7 

.5 

.2 

.5 

.5 

.1 

.3 

.2 

.2 

1.0 

1.0 

.1 

.4 

3.0 


.5 

.7 

1.2 

1.0 

1.0 

.5 

.1 

.5 

.8 

1.0 

.5 

1.0 

.2 

.1 

1.0 


.5 

1.0 

1.0 

.1 

.2 

.5 

.5 

.5 

.2 

.1 

.2 

.5 

2.0 

1.0 

.4 


.5 

1.0 

.6 

1.0 

1.0 

2.0 

.5 

.2 

1.0 

3.0 

.2 

.5 

.5 

.5 


1.0 

.5 

.4 

.9 

.5 

.5 

1.0 

.5 

.2 

.9 

.1 

.5 

.5 

1.5 


.5 

1.0 

.1 

.5 

1.5 

.5 

1.0 

1.0 

.Z 

.5 

.9 

.2 

.5 

1.0 


.2 

.5 

.5 

.1 

1.0 

.2 

.5 

1.0 

1.0 

.2 

1.2 

.1 

.2 

1.0 
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' TABLE  II 

None  Grouped  Frequency  Distribution  of  Repair  Times 


Class  Intervals 

Frequency 

Cumulative  Frequency 

0 - .1 

26 

26 

.1  - .2 

54 

80 

.2  - .3 

13 

93 

.3  - .4 

8 

101 

.4  - .5 

249 

.5  - .6 

9 

258 

.6  - .7 

5 

263 

.7  - .8 

2 

265 

.8  - .9 

4 

269 

.9  - 1.0 

67 

336 

• 

1 

O 

0 

336 

1.1  - 1.2 

2 

338 

1.2  - 1.3 

1 

339 

1.3  - 1.4 

0 

339 

1.4  - 1.5 

8 

347 

1.5  - 1.6 

0 

347 

1.6  - 1.7 

0 

347 

00 

• 

1 

e 

0 

347 

1.8  - 1.9 

0 

347 

1.9  - 2.0 

9 

356 
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TABLE  II  (Continued) 


None  Grouped  Frequency  Distribution  of  Repair  Times 
Class  Intervals  Frequency  Cumulative  Frequency 


2.0  - 2.1  - . 

0 

356 

2.1  - 2.2 

0 

356 

2.2  - 2.3 

0 

356 

2.3  - 2.4 

0 

356 

2.4  - 2.5 

2 

358 

2.5  - 2.6 

0 

358 

2.6  - 2.7 

0 

358 

2.7  - 2.8 

0 

358 

2.8  - 2.9 

0 

358 

2.9  - 3.0 

3 

361 

3.0  - 3.1 

0 

361 

3.1  - 3.2 

0 

361 

3.2  - 3.3 

0 

361 

3.3  - 3.4 

0 

361 

3.4  - 3.5 

1 

362 

3.5  - 3.6 

0 

362 

3.6  - 3.7 

0 

362 

3.7  - 3.8 

0 

362 

3.8  - 3.9 

0 

362 

3.9  - 4.05 

4 

366 
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MrEW«L  SlilDS  F^ESJEVCr  GR4P-I 


EACH  X REPRESENTS  2.A66667  OBSERVATIONS  Figure 


K 

K 

*c  ye 
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M K 
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K W 
K K K 
K M 
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MPlINT.  REPRIR  TIMES  (HRS) 

Each  5th  Observation  is  Represented  by  + 


TABLE  III 


I 


-Gsoupad  Frequency  Distribution  of  Repair  Times 


Cl«iss  Intervals 

Frequency 

Cumulative  Frequency 

0 - 

.15 

26 

26 

.15  - 

.25 

54 

80 

.25  - 

.75 

183 

263 

.75  - 

1.25 

75 

338 

1.25  - 

1.75 

9 

347 

1.75  - 

2.25 

9 

356 

2.25  - 

4.05 

10 

366 

EXPONENTIAL 


[ 

i 


! interval 

IHRS) 

OBSERVED 

EXPECTED 

CO-E  I^^E/E 

SUN 

.00 

- 

.15 

26.00 

74.07 

31.19 

31.19 

.15 

« ’ 

- 

.25 

54.03 

40.85' 

4.23 

35.43 

.25 

- 

.35 

13.33 

35.13 

13.94 

49.37 

.35 

- 

.45 

8.00 

30.22 

16.34 

65.71 

.45 

- 

.55 

148.00 

25.99 

572.80 

638.51 

! -55 

- 

.65 

9.33 

22.35 

7.98 

646.49 

1 .65 

- 

.75 

5.00 

19.22 

10.53 

657.01 

.75 

- 

.35 

2.30 

16.53 

12.78 

669.79 

.85 

1 

- 

.95 

4.00 

14.22 

7.35 

677.13 

.95 

- 

1.35 

67.03 

12.23 

245.24 

922.38 

1.35 

- 

1.15 

.30 

10.52 

10.52 

932.90 

1.15 

- 

1.25 

2.00 

9.05 

5.49 

938.39 

1.25 

- 

1.35 

1.00 

7.78 

5.91 

944.30 

1.35 

- 

1.45 

.33 

6.69 

6.69 

950.99 

l.<i5 

- 

1.55 

8.00 

5.76 

.87 

951.87 

1.55 

- 

1.75 

.30 

9.21 

9.21 

961.07 

1.75 

- 

1.95 

.00 

6.81 

6.81 

967.89 

1.95 

- 

2.15 

9.30 

5.04 

3.11 

971.00 

2.15 

- 

2.45 

.33 

5.21 

5.21 

976.21 

2.45 

- 

4.35 

10.00 

8.29 

.35 

976.56 

Table  IV 
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EXPONENTIAL 


INTERVAL 

IHRSI 

OBSERVED 

EXPECTED 

C0-E>«*2/E 

SUN 

.00 

- 

.15 

26.00 

76.07 

31.19 

31.19 

.15 

- 

.25 

56.00 

60.85 

6.23 

35.63 

.25 

- 

.75 

1B3.00 

132.92 

18.87 

56.30 

.75 

- 

1.25 

75.30 

62.55 

2.68 

56.77 

1.25 

- 

1.75 

9.00 

29.66 

16.19 

70.97 

1.75 

- 

2.25 

9.00 

13.86 

1.70 

72.67 

2.25 

- 

^.35 

10.30 

11.50 

.20 

72.86 

Table  V 
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WEIBULL 


INTERVAL 

(HRS) 

OBSERVED 

EXPECTED 

lO-E l••2/E 

SUN 

o 

o 

• 

- 

.15 

26.00 

61.47 

20.47 

20.47 

.15 

- 

.25 

54.03 

39.86 

5.02 

25.48 

.25 

- 

.35 

13.00 

36.13 

14.81 

40.29 

.35 

- 

.'«5 

8.00 

32.13 

18.12 

58.41 

.45 

- 

.55 

148.00 

28.25 

507.58 

565.99 

.55 

- 

.65 

9.00 

24.64 

9.93 

575.92 

.65 

- 

.75 

5.00 

21.36 

12.53 

588.45 

.75 

- 

.95 

2.00 

18.43 

14.65 

603.10 

.85 

- 

.95 

4.00 

15.84 

8.85 

611.95 

.95 

- 

1.15 

67.00 

13.56 

210.53 

822.48 

1.05 

- 

1.15 

.00 

11.58 

11.58 

834.06 

1.15 

- 

1.25 

2.00 

9.86 

6.27 

840.33 

1.25 

- 

1.35 

1.00 

8.38 

6.50 

846.83 

1.35 

- 

1.45 

.00 

7.10 

7.10 

853.93 

1.45 

- 

1.55 

8.00 

6.01 

.66 

854.59 

1.55 

- 

1.65 

.00 

5.08 

5.08 

859.67 

1.65 

- 

1.35 

.00 

7.89 

7.89 

867.56 

1.85 

- 

2.)5 

9.03 

5.58 

2.10 

869.65 

2.05 

- 

2.35 

.00 

5.42 

5.42 

875.07 

2.35 

- 

4.35 

10.00 

7.13 

1.16 

876.23 

I 

TaU.e  VI 


1 


INTERVAL 

IHRS) 

aSSERVEO 

EXPECTED 

lO-E  l••2/E 

SUN 

.0 

- 

.15 

26.33 

61.67 

20.67 

20.67 

.15 

- 

.25 

56.33 

39.86 

5.02 

25.68 

.25 

- 

.75 

183.00 

162-52 

11.50 

36.98 

.75 

- 

1.25 

75.00 

69.26 

.67 

37.66 

1.25 

- 

1.75 

9.33 

30.86 

15.68 

52.96 

1.75 

- 

2.25 

9.00 

13.11 

1.29 

56.23 

2.25 

- 

A. 35 

10.00 

8.62 

.22 

56.65 

39 


GAMMA 


INTERVAL 

(HRSi 

OBSERVED 

EXPECTED 

(0-E l••^/E 

SUM 

• 'J  0 

- 

• 1 5 

26.00 

12.9A 

13.17 

13.17 

.15 

- 

.25 

5A.00 

28.3A 

23.23 

36.  AO 

.25 

- 

.35 

13.00 

39.29 

17.59 

53.99 

.35 

- 

.«5 

8.00 

A3.7A 

29.20 

83.19 

.^5 

- 

.55 

148.00 

A3.1A 

25A.88 

338.08 

.55 

- • 

.55 

9.00 

39. A 1 

23. A7 

361. 5A 

.65 

- 

.75 

5.00 

3A.15 

2A.88 

38  6.  A3 

.75 

- 

.35 

2.00 

28.A7 

2A.61 

All.OA 

.85 

- 

.75 

A.  00 

23.0A 

15. 7A 

A26.77 

.9  5 

- 

1.  '5 

67.00 

18.22 

130.63 

557.A0 

1.05 

- 

1 .15 

.00 

1A.13 

1A.13 

571. 5A 

1.15 

- 

1.2  5 

2.00 

10.80 

7.17 

578.70 

1.25 

- 

1.35 

1.00 

8.1A 

6.26 

5BA.97 

1.35 

- 

1.45 

.00 

6.07 

6.07 

591. OA 

1.45 

- 

1*55 

8.00 

7.77 

.01 

591.05 

1.65 

- 

A.  j5 

19.00 

8.35 

13.58 

60A.63 

TaTal'e  VIII 
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*■*'  - ■ .. 


V • 


GAMMA 


INTERVAL 

(HRS) 

OBSERVED 

EXPECTED 

lO-E  l••^/E 

SUN 

.00 

- 

• 13 

26.00 

12.94 

13.17 

13.17 

.15 

- 

.25 

54.00 

28.34 

23.23 

36.40 

.25 

- 

.73 

IS3.00 

199.73 

1.40 

37.80 

.75 

- 

1.25 

75.00 

94.66 

4.08 

41.89 

1.25 

- 

1.75 

9.00 

24.37 

9.69 

51.58 

1.75 

- 

4.  >3 

19.00 

5.96 

28.51 

80.08 

Tame  IX 
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S OCT  7o  PAGE 


TRUNCATED  NORMAL 


IIiTERVAL 

(I-RS) 

D-OSf  AVEO 

EXPECTEO 

(0-E)**2/E 

SUM 

.00 

- 

.15 

26.00 

67,44 

25.46 

25,46 

.’l5 

- 

.25 

54.00 

38,31 

5.94 

31.41 

,25 

- 

.35 

13.00 

34,40 

13.31 

44,72 

.35 

- 

.45 

8.00 

30.39 

16,50 

61.22 

.^5 

- 

.55 

148.00 

26,77 

549,03 

610.25 

.55 

m 

.65 

9.00 

23,50 

3,95 

619.20 

.65 

- 

.75 

5.00 

20,57 

11.79 

630.98 

.75 

- 

.35 

2.00 

17,95 

14.17 

645,15 

.35 

m 

.95 

4.00 

15.61 

3,64 

653.79 

.95 

- 

1.05 

67.00 

13,54 

211.15 

364,94 

1.05 

- 

1.15 

.00 

11,70 

11.70 

876 • 64 

1.15 

- 

1.25 

2.00 

10.08 

6.48 

983.12 

1.25 

- 

1.33 

1.00 

3.66 

6,78 

389.90 

1,35 

- 

1.45 

.00 

7.42 

7.42 

897.31 

l.<*b 

- 

1.55 

8.00 

6,33 

.44 

897,75 

1.55 

- 

1.65 

.00 

5.39 

5.39 

903,14 

1.05 

«• 

1.75 

.00 

4,57 

4.57 

907.71 

1.75 

1.95 

.00 

7.12 

7.12 

914.83 

1.95 

- 

4.05 

19.00 

15,60 

.74 

915.57 

TABLE  X 
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8 OCT  fh  ?>.,or; 


TRUNCATED  NORMAL 


INT 

: 3.  V .t  L 

(ii7S  ) 

a-SE»vEo 

Expected 

(0-E)»>=<'2/E 

SUf' 

9 - J 

- 

.15 

26.  on 

(•  7.  -.4 

25.46 

25.46 

.It)’ 

.25 

JO 

33. ai 

5,04 

31.41 

.23 

- 

.75 

I '3.00 

135.64 

16,54 

47,94 

.75 

- 

1.25 

75.00 

, 36 

.54 

48,49 

1.23 

- 

1.75 

9,u0 

32.36 

16.36 

65.35 

1.73 

- 

1.95 

9.00 

7,12 

,50 

65.85 

l."^3 

- 

•>.J5 

10.00 

15, 6J 

2.01 

67,86 

TABLE  XI 
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LOGNORMAL 


INTER)/ AL 

(HRS) 

OBSERVED 

EXPECTED 

IO-£»**2/E 

SUM 

. .0 

- 

.15 

26.00 

22.21 

.65 

.65 

.15 

- 

.25 

54.00 

46.53 

1.20 

1.85 

.25 

- 

.35 

13.33 

50.81 

28.13 

29.98 

.35 

- 

.45 

8.30 

45.53 

30.94 

60.92 

.45 

- 

.35 

148.33 

37.92 

319.49 

380.41 

.55 

- 

.55 

9.00 

30.66 

15.30 

395.72 

.65 

- 

.75 

5.30 

24.51 

15.53 

411.25 

.75 

- 

.35 

2.33 

19.55 

15.75 

427.00 

.85 

- 

.95 

4.00 

15.61 

8.64 

435.64 

.95 

- 

1.35 

57.00 

12.52 

237.19 

672.82 

1.05 

- 

1.15 

.33 

10.08 

10.08 

682.90 

1.15 

- 

1.25 

2.30 

8.17 

4.66 

687.56 

1.25 

- 

1.35 

1.30 

6.65 

4.80 

692.37 

1.35 

- 

1.45 

.00 

5.45 

5.45 

697.82 

1.45 

• 

1.55 

8.00 

8.21 

.01 

697.82 

1.65 

- 

1.35 

.33 

5.68 

5.58 

703.51 

1.85 

- 

2.15 

9.00 

5.56 

2.13 

705.63 

2.15 

4.35 

10.00 

9.18 

.07 

705.71 

•.Tame  ni 
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LOGNORMAL 


INTERV  AL 

(HAS) 

OBSERVED 

EXPECTED 

I0-EI**2/E 

SUN 

, 0 

- 

.15 

2b. OQ 

22.21 

.65 

.65 

.15 

- 

.25 

Sb.OO 

46.53 

1.20 

1.85 

.25 

- 

.75 

183.00 

189.44 

.22 

2.07 

.75 

- 

1.25 

75.30 

65.92 

1.25 

3.32 

1.25 

- 

1.75 

9.00 

23.41 

8.87 

12.18 

1.75 

- 

2.25 

9.00 

9.47 

.02 

12.21 

2.25 

- 

4.15 

10.00 

7.86 

.58 

12.79 
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Commander 
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ATTN:  DRSAV-D 
PO  Box  209 

St.  Louis,  MO  63166  1 
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